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Abstract 

We compute the two loop coefficient in the relation between 
the lattice bare coupling and the running coupling defined through 
the Schrodinger functional for the case of pure SU(2) gauge the- 
ory. This result is needed as one computational component to 
relate the latter to the MS-coupling, and it allows us to imple- 
ment 0(a) improvement of the Schrodinger functional to two- 
loop order. In addition, the two-loop /3-function is verified in a 
perturbative computation on the lattice, and the behavior of an 
improved bare coupling is investigated beyond one loop. 



1 Introduction 



The lattice formulation of pure gauge theories and QCD enables one, in 
principle, to connect the non-perturbative (low energy) regime to the per- 
turbative (high energy) regime. An important example is the computation 
of the running coupling at short distances in some scheme, say minimal sub- 
traction (MS), with the scale set by a low energy quantity of the theory. 
At the outset there are two problems that make this computation practi- 
cally difficult. One is the apparent need for large lattices to span a wide 
range of length scales. This problem can be circumvented by the use of 
finite-size scaling as was first demonstrated for non-linear sigma models in 
two-dimensions [1]. The second problem is to find a physical renormalized 
quantity to serve as a running coupling, that has a number of favorable prop- 
erties [2]. Firstly, it should allow for an accurate numerical estimation on 
a range of lattices. Secondly, it should exhibit small cut-off dependence, so 
that a reliable and precise extrapolation to the universal continuum values 
is possible from medium size lattices. Thirdly, although the definition has to 
be non-perturbative, the quantity must have a manageable perturbation ex- 
pansion at high energy. This is required to perturbatively relate the coupling 
chosen for the simulation to other more conventional ones at high energy. 

In refs. [3], [4], [5] the coupling g 2 (L), based on the Schrodinger functional 
with an induced abelian background field, was found to fulfill these require- 
ments to a sufficient degree. In the present paper, for the case of SU(2), we 
work out the 2-loop contribution to g 2 on the lattice in terms of the bare 
coupling (/q = 4//3. Recently [6] a calculation of the relation between and 
^llg was reported and combined with our result. We have thus succeeded 
in connecting g 2 and g^, at our highest energy with an accuracy needed to 
match the improved statistical precision on parallel computers [7]. In addi- 
tion, the present computation alone allows us to implement 0(a) Symanzik 
improvement of our action and coupling to 2-loop order. As a by-product, 
we investigate the expansion of g 2 in terms of a simple plaquette-improved 
bare coupling [8] to one more order than before. This is highly welcome in- 
formation, as up to now, only the smallness of 1-loop contributions in terms 
of such couplings could be verified for some quantities. 

The organization of the paper is as follows. In sect. 2 we briefly remind 
the reader of the definition of g 2 in terms of the Schrodinger functional and 
set up the necessary expansions. In sect. 3 we list perturbative results for our 
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coupling on a series of lattices and analyze their cut-off dependence. Up to 
here the paper should be readable in a self-contained fashion for a reader, who 
just wants to understand the problem and the result. The remainder contains 
more technical material of interest to those who plan a similar computation. 
In that part we assume that the reader has a copy of [3] at hand, and we 
refer to individual equations and sections of this reference. In sections 4 and 
5 we discuss the structure of the gauge fixed action and the construction 
of propagators in the background field. In sect. 6 all terms required at 2- 
loop order are collected and the optimization of the computational speed is 
discussed. 



2 Running coupling and Schrodinger func- 
tional 

In this section we summarize the steps in defining the Schrodinger functional 
based running coupling g 2 proposed in [3] and define its perturbative expan- 
sion. We use lattice units, setting the spacing a to unity, and choose a box 
with temporal and spatial extent L. 

The Schrodinger functional is the kernel for euclidean propagation from 
some field configuration at time x° = to some other configuration at x° = L. 
In the path integral formalism, this becomes the partition function on an L 4 
box with fixed boundary conditions at x° = 0, L and periodic ones in the 
three space directions. 

A lattice gauge field U is an assignment of a link variable U(x, fj,) £ SU(2) 
to every pair (x, x + ft) of lattice points. In the path integral to be defined, 
all gauge fields in the interior fluctuate freely, and the ones on the spatial 
links in the time slices x° = 0, L remain fixed at 

U(x,k)\ x o =0 = exp{C fc }; U(x,k)\ x o =L = exp{C' k }; k = 1,2,3, (2.1) 

with 

C k = ~r , C^—j-r, r = ^ _J. (2.2) 

These boundary conditions contain one dimensionless free parameter with 
range [3] < r\ < ir. The Schrodinger functional is now given by 

Z( V ,L) = ( D[U]e- s V\ D[U] = Y[dU(x,fj.), (2.3) 
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where one integrates over all interior lattice gauge fields with the fixed bound- 
ary values given by eq. (2.1). The action S[U] is taken to be 

S l U ] = \J2 «v(s)tr{l " U^(x)}, (2.4) 

where the sum runs over all oriented plaquettes on the lattice and U fJil/ (x) 
denotes the parallel transporter around the plaquette. The weight w fJil/ (x) is 
given by 

{\ c s{go) f° r spatial plaquettes at x° = or x° = L, 
c t(go) f° r time-like plaquettes touching the boundary, (2-5) 
1 elsewhere. 

In the classical limit the choice c s (0) = c t (0) = 1 eliminates 0(a) lattice 
artifacts. Beyond it, the functions c s (go) and c t (go) are chosen to minimize 
cut-off effects. For the abelian boundary fields (2.2), c s (go) is irrelevant and 
will hence not be discussed further (see however [3] for 1-loop values of both 
weights). For perturbative Symanzik improvement we insert 

c t (g ) = 1 + c^gl + c[ 2) g* + ■■■ (2.6) 

(2) 

and will report the 2-loop coeffient c\ in the next section. 
With the present boundary conditions the gauge field 

V(x,0) = l; V(x,k) = exp{[x C' k + (L-x )C k ] /L} (2.7) 

and its gauge transforms minimize the classical action with a value 

T ( V ,L) = 9 2 S[V]\ g2o=o = 24L 4 sin 2 {^^} • (2.8) 

The response to an infinitesimal change in r\ defines the renormalized coupling 
proposed in [3], 



=2,M r 0(^) 



(2.9) 

?7=7r/4 



V'( V: L 
where 

T(ri,L) = -In Z( V ,L) (2.10) 
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is the effective action. The primes denote derivatives with respect to r\. Here, 
and in all simulations sofar, the value r\ = ir/4 has been chosen [3], which 
gives a good signal in the Monte Carlo simulation. 

In perturbation theory in go the expansion of the effective action 

r = ^ - 2 r + r 1 +^ 2 r 2 + --- (2.11) 

yields 

g 2 (L) = gl + mi {L)g* + m 2 (L)g 6 + • • • , (2.12) 

where 

p/ p/ 

mi = --f; m 2 = --f + m 2 . (2.13) 
1 o 1 o 

The computation of rri2(L) is the subject of this article. 

3 Results 

The perturbative expansion of the Schrodinger functional on the lattice with 
its finite induced background field can only be evaluated numerically. We 
carry this out for a range of lattices up to L = 32 and analyze the asymp- 
totic L-dependence of the various expansion coefficients of interest. The 
general technique follows closely what is described in [10], and details of our 
calculation are given in the next sections. Here we only mention that the 
numbers we report here underwent a battery of tests like independence of 
a free gauge fixing parameter and verification of symmetries before actually 
using them to shorten the sums. Moreover, the two authors created two 
completely independent programs (separated by the atlantic ocean), which 
yielded identical results up to L = 8. Beyond this size, we only ran the 
program that has been optimized for speed. The total effort is on the order 
of ~ 1000 h on a Sun 10-41. 

3.1 Coefficients for the coupling from L = 4 . . .32 

To exhibit the dependence on improvement we parameterize rai(L), rri2(L), 
mi = itt-i + c^m\ (3-1) 

a , (1) b , r (1)1 2 c , (2) d , 2 r,\ 
1712 = nT-2 + C t m 2 + \ C t m 2 + C t m 2 + m l (3-2) 
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where the occurring powers of cj 1 ' follow from a glance at the structure of 
the perturbative terms (sect. 6). Some of the contributions are simple and 
can be given in closed form 

2 

m\ = m\ = — — (3-3) 



L 

1 + cos" 2 



7r — 2r/ 



(3.4) 



Numerical results for the remaining terms are listed in table 1. Values for 
toi(L) were partly quoted 1 in [3] and are extended here for completeness. 
They are known in extended precision (128 bits), and all digits in the table 
are significant. The 2-loop computation was done in double precision (64 
bits), and the truncation is such that the expected error is in the last digit. 

3.2 Perturbation theory of the step scaling function 

In our previous work [1] - [5] we introduced the step scaling function to 
characterize the evolution of a renormalized coupling under a finite rescaling 
factor s. For g 2 it is defined in the continuum by 

a(s,u) = f(sL) \ u= - 2[L) . (3.5) 

On the lattice, this universal function [3] emerges in the limit from pairs of 
lattices of sizes L, sL which diverge in lattice units, 

ais.u) = lim T,(s,u,l/L). (3.6) 

l/L-s-0 

Both S and a possess a perturbative expansion in u. Of particular interest for 
the numerical work is the perturbative expansion of the effect of the lattice, 

S{u, 1/L) = S(W/^)-g(2,tO = Si{Uj 1/L) u + H ^ 1/L) u2 + _ _ (3 7) 

The coefficients ^1,^2 can be straightforwardly derived from mi, 1712 and are 
collected in table 2. We notice that the 1-loop improved and the 2-loop 
improved 8^ decay roughly proportional to 1/L 2 as expected. 



""in [3], mi(L) was defined not to include the contribution from 1-loop improvement 
and corresponds to mJ(L) here. 
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U 


m a 

1 1 0-1 




m b 
,n 2 


A 
L ± 


31 301 991 1 fi^R31 8 


D1 84Q43Q1 1 40fi97fi 


0fifi4fi8304071 1 000 


K 
o 


335655931 7389083 


0901 599604839990 

.\JLi\J ±OLi*y\J\J^:OOLiLiLi\J 


0630548439988880 


u 


3^14990301 8055683 

,00^:LiLi\JO\J ±0\JOO\JOO 


D91 4004669936067 


0504509734930983 

.\JOZ>^:OZ> Li t O^LiOVJLiOO 


7 


37DD1 769371 831 Q5 

.Of UU1 I UZiO 1 IOOIju 


0993569763764735 


0560994083970087 

.UJUU LiLi'-iVJOO Li t Z>\JO t 


8 

o 


383704988961 791 9 

.OOOI U^:^OO^Ul l Li X Li 


0931 1 41 31 397989 


059839790059609 

.\JO LiOO Li 1 LiWOLi\j\JLi 


Q 


30^17330981 47R000 


09373R0Q7R4M 


04001 QK.QC11 9^1409 


1 n 


40R441 R740fi03030 


0949R1 9M R479RR 


047970^139908080 


1 1 


41 6D77775D7D599Q 

.t:1UUI l l t OVJ t VJO LiLiZ? 


D9471 4R38901 379 


044809009849661 


1 9 

± Li 


494831 D66Q1 1 3876 


0951 1 9969750696 


049734965040504 


1 3 


4398461 8Q54D334D 


0954678701 1 3941 


0407776681 80908 

.Ut:UI l t \J\JO ±Oc> LiVJO 


1 4 


44093*i77Q41 9R31 R 


09^178771 83934R4 

.UZiO /Oil IOOZiOtiUt: 


03RQQRQ*iQ971 1 *i3 


1 ^ 


4470RRQR7fi3Q7QR7 


09807873^070^148 


03737RR09004397 


1 fi 


4^3477440^677^1 fi 


09634565788463 


03^8Q1 783Q69939 


1 7 


4^1045051 04479748 


0965091 631 3431 


03459893365001 6 

.\JO^OLiOLiOO\JO\JZ>±\J 


1 8 

J. O 


465D83391 69Q1 1 34 


096891 1 691 9734 


033971 039739646 


1 Q 

X C* 


47D388Q8591 41 854 

.Til \JOOOc>00 Li IIIOiJt: 


09703400336057 


0391 1 081 1 438987 


90 


47^41 030RR1 01 1 9^ 


09793^61 R4R39 

,\j Li t LiOOO\J X OttUOZi 


031 0344R4031 R1 4 

.U01U0t:t:Ut:U01 U It: 


91 

z< X 


48D1 76DD63D68773 


09749443844983 


0300338074831 70 


LiLi 


48471 n6884549DD1 


097609039751 90 

.\JLi 1 \J\J LiO O Li 1 O j-LiZJ 


0901 01 31 5688398 


93 

LiO 


48QD355885758341 


0977791 4790876 

.\JLi t t t Li Xt: i LiuO t U 


098990036049600 

,\J Lj\jLjLjVVO\JV^iLi\Jij\J 


24 


.4931691576726372 


.02793301855603 


.027413840712453 


25 


.4971275264639118 


.02808634837864 


.026647839420229 


26 


.5009248770736207 


.02823282841524 


.025927371235302 


27 


.5045737434715104 


.02837305927392 


.025248410342669 


28 


.5080852563831408 


.02850756573057 


.024607391692524 


29 


.5114693444997410 


.02863680891327 


.024001147727234 


30 


.5147349009901648 


.02876119615134 


.023426855017935 


31 


.5178899222329467 


.02888108900318 


.022881989110948 


32 


.5209416241340969 


.0289968098474 


.022364286191250 



Table 1: List of L-dependent coefficients m",™^ and m h 2 . 
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L 


4 0) 


s[ 1] 


°2 


°2 


°2 


4 





00631 


-0 


00726 





00255 





00224 


-0 


00064 


5 





00641 


-0 


00445 





00204 





00194 


-0 


00036 


6 





00623 


-0 


00282 





00163 





00168 


-0 


00024 


7 





00584 


-0 


00192 





00132 





00147 


-0 


00017 


8 





00540 


-0 


00139 





00109 





00131 


-0 


00012 


9 





00497 


-0 


00106 





00091 





00118 


-0 


00009 


10 





00459 


-0 


00084 





00077 





00108 


-0 


00007 


11 





00425 


-0 


00068 





00066 





00099 


-0 


00006 


12 





00396 


-0 


00056 





00057 





00091 


-0 


00005 


13 





00370 


-0 


00048 





00050 





00085 


-0 


00004 


14 





00347 


-0 


00041 





00044 





00079 


-0 


00003 


15 





00327 


-0 


00035 





00038 





00074 


-0 


00003 


16 





00309 


-0 


00031 





00034 





00070 


-0 


00002 



Table 2: Perturbative lattice effects in the step scaling function, eq. (3.7). 
The superscripts in brackets here refer to 0, 1 and 2-loop 0(a) improvement. 

3.3 Asymptotic analysis 

In the spirit of ref. [10] we now derive the asymptotic expansion of mi and m,2 
in terms of logarithmically modified powers of the lattice spacing, L~ m (ln L) n , 
with n < 2 for our 2-loop calculation. As stressed in [10], the generally 
best way to access the successive terms in the expansion is to recursively 
take discrete L-differences in such a way, that the desired term survives, 
while all the more leading ones and some of the next subleading ones are 
canceled. Over a range of L the approximate constancy (up to uncancelled 
subleading terms) of this combination has to be checked. In this process 
both the available range of L as well as the number of significant digits 
shrink, which limits the degree of cancellation possible. An analysis of both 
systematic [10] and roundoff errors has been implemented. It also has to be 
born in mind, that the expansion in 1/L is expected to be asymptotic only. 
We saw indications that for L < 32 the degree in 1/L that we could access 
with our 64-bit precision, roughly matches with the number of terms where 
the expansion makes sense. This means, that to reduce the errors in the 
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numbers to follow, we would need both more digits and larger L. Neither 
of them alone would help significantly. An extension in this direction — 
fortunately not urgently needed — would seem rather difficult, as there seem 
to be hardly any compute-resources available that are efficient for large scale 
number crunching in 128-bit precision. 
The result for m\ in [3] is 

m%(L) = 2b In L + 0.202349(3) - 0.1084(11)-^ + (^Jjrj > ( 3 - 8 ) 
where the first of the universal coefficients of the /3-function 



60 = 24^' 61 = 9^' (3 - 9) 



enters. Demanding the absence of a 1/L artifact in mi fixes the 1-loop 
improvement coefficient [3] 



-0.0543(5). (3.10) 



For we find 



In L ( 1 (In L 



\2 N 



m a JL) = 2&i In 1 + 0.016069(3) + 0.011(1)— + -, v T ' . (3.11) 

L \L L 2 J 

The absence of two terms is noteworthy here. If we write down the Callan- 
Symanzik equation for V and insert eq. (2.11) and the perturbative f3- 
function, it is trivial to infer that m%(L) should have no (In L) 2 divergence, 
as we indeed find. In fact, we observe the leading 2&i In L behavior with a rel- 
ative accuracy of 10~ 3 . To our knowledge, this is actually the first time that 
the 2-loop coefficient of the /3-function is observed in lattice perturbation 
theory. The exact logarithm is then subtracted for the following analysis. 
The absence of the most leading artifact term oc (In L) 2 /L is expected from 
tree-level improvement, and we can bound its coefficient by ~ 10~ 3 . The 
coefficient of 1/L could not be determined reliably from m^L) alone. 

In a next step, we analyze m h 2 (L), which is composed of 1-loop diagrams, 
and determine the leading divergence 

c{ 1} m b 2 (L) = -0.01008(5)^ + (j) ■ ( 3 - 12 ) 
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We here used an L-dependent estimate for c\ from a suitable combination of 
the 1-loop results for mi, and we analyzed the L-dependence of the product as 
a whole. We see that the cancellation of In L/ L with the term in (3.11), which 
is necessary for full 0(1/ L) improvement, is consistent with the data, with 
the more precise value coming from m h 2 (L). Amusingly, this also constitutes 
another check on cj 1 '. 

We finally analyze the combination 

m a 2 + c{ 1} m b 2 +\c{ 1} ] m c 2 -2b 1 lnL = 0.016069(3) -0.0230(9)- + 

(3-13) 

All terms of first order in the lattice spacing are now canceled in m,2 by 
choosing 

c[ 2) = -0.0115(5) (3.14) 
as the 2-loop improvement coefficient. 

3.4 Two-loop relations between couplings 

In this subsection we discuss the perturbative relation between our coupling 
ct(q) = g 2 (L)/Air with q = 1/ L in the continuum and the bare coupling on the 
lattice 2 Q!o = g^/Air. This relation is expected to be free of large logarithms 
if q is of order a -1 , where a is the lattice cut-off length 3 associated with go- 
lf is obtained by neglecting all inverse powers of L and then setting L = a, 
for instance, in the formulas of the previous subsection, 

c^a" 1 ) = a + 2.543 a 2 + 9.00 a 3 + O(ag). (3.15) 

We see rather large coefficients here confirming the common experience, that 
cto is not a good expansion parameters for physical quantities. For f3 = 2.85 
we have cto ~ 0.11 and hence corrections of 28% and 11%. It has been 
heuristically argued [8], and also observed at the 1-loop level, that a better 
coupling is given by 

a = a /P, (3.16) 

2 These results were also summarized in [9]. 
3 Only for this subsection do we reintroduce a. 
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where P is the mean plaquette at the corresponding bare coupling (&q ~ 0.16 
at (3 = 2.85). With its expansion [11], [6] 

P = 1 ~TQ g2 °~ °- 00946141 4o + OG/o 6 ), (3-17) 

(3.15) resums to 

a (a _1 ) = a + 0.1866 a 2 + 1.08 &l + 0(a£). (3.18) 

The 1- and 2-loop corrections are about 3% each. In view of this, the small- 
ness of the 1-loop term looks somewhat accidental. 

We realize that the choice of exactly 1/a as the scale for a instead of 
some other number of the same order, like for instance 2ir / a suggested by 
momentum discretization, is arbitrary. This clearly limits the significance 
of the size of the 1-loop term. With the 2-loop term at hand we have the 
possibility to fix this freedom in a more compelling way. Choosing the scale 
to minimize the 2-loop coefficient we derive 

a(1.45a _1 ) = a - 0.2460ag + 0.891(l)ajj + O(a2), (3.19) 
a^O^la" 1 ) = a - 0.2460ag + 1.226(l)ag + O(ag). (3.20) 

The size of the coefficients is similar in both cases here. In fact, due to 
the larger value of &q : this scheme has bigger corrections at f3 = 2.85. The 
identical 1-loop coefficients at the minimizing scale are no coincidence. They 
are given analytically as —2irbi/bo. Another possibility is to choose a scale 
where the 1-loop correction vanishes. This results in the qualitatively similar 
relations 

a(1.17a _1 ) = a + 0.951(l)ag + 0{6t%), (3.21) 
a(8.83a _1 ) = a + 1.287(l)ajj + O(c^). (3.22) 

In the light of these 2-loop results there arises the possibility that the 
main difference between the standard bare coupling and the modified bare 
coupling &q consists of the scale that it is best used for. This point has been 
further clarified by comparing with our new high precision data in ref. [7]. 

The relation between a and cty^, is discussed in [6]. The extension of the 
present calculation to SU(3) is planned. 
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4 Expanded gauge fixed action 



In the remaining sections we report details of our computation. It will be of 
particular interest to those readers embarking on a similar lattice perturba- 
tive calculation. We shall refer from now on to section and equation numbers 
of ref. [3]. 

The perturbative expansion on the lattice requires to fix the gauge. This 
amounts to supplementing the pure gauge action S in (2.4) by a gauge fixing 
term S g { together with a Fadeev-Popov ghost term Spp involving additional 
fermionic fields. We shall now discuss these contributions as series in go up 

to oG/o 2 ). 

4.1 Some notation 

We specialize the discussion in [3] to SU(2). In its Lie algebra we use the 
basis T a , a = +, — , 0, 

T± = ^(nisra) (4.1) 
Zi 

T° = Ir 3 (4.2) 

with the Pauli matrices r^. For vector potentials q^x) = q tl>(T ( y x)T a and 
r ni x ) = r fj,,a{x)T a (sum over a) the scalar product is then given as 

(l, r ) = J2i 2< in,+ ( x ) r ^-i x ) + 2 %A X ) r nA x ) + %fi{ x ) *V,o(a5)}, (4-3) 

and there is a corresponding formula for scalar fields without the //-sum. We 
define the raising of indices a by = 2q^^, q® = q^fl an d may then write 

(9,r)= £ lZ( x )r,A x )- (4-4) 

4.2 Fields and boundary conditions 

The abelian background field V in (2.7) is parameterized by 

V(x,k) = e-^ T ° (4.5) 
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with 

(3{x ) = -{2lL 2 )[vL + (tt - 277)30]. (4-6) 
Its constant field strength is 

= e^ T ° (4.7) 

with 

7 ^ = _ 7 ^ = ((2/ J C 2 )(7r-2r / ) for p = 0, u = 1,2,3 (4 g) 
I else 
It is convenient to parameterize the fluctuations of the lattice gauge field 
around V by 

U{x,n) = eyLv{g q^{x)}V{x, fj.) (4.9) 

with in the Lie Algebra of SU(2). Our gauge fixing will lead to all modes 
of being quadratically damped. Hence in perturbation theory we only deal 
with infinitesimal goq^. 

In the Schrodinger functional we have to integrate over qo t(T (x) for xq = 
. . . L — 1 and over qk,<r{%) for xq = 1 ... L — 1 {k = 1, 2, 3). As discussed 
in detail in [3], sect. 6, the notation is considerably simplified, if we define 
additional fields by boundary conditions. Since all fields are periodic in space, 
we can discuss them in terms of the Fourier transforms ^^(xq, p), 

<Mo,p) = <M£,p) = o (4.10) 

qoAL,P) = <?cv(£-l,p) (4.11) 
qoA-hp) = <?cv(0,p) if (a,p)^(0,0) (4.12) 
<?o,o(-l,0) = 0. (4.13) 

The Haar measure in the q^ parameterization results in a local contribution 
to the action S m , 

J] dU(x, (j.) = n dq^{x) exp{-S m } (4.14) 

with 

2 



^ = ^E?,V + o(4 (4-15) 

The Fadeev-Popov determinant is represented with Lie algebra valued 
ghost fields c(x) and c(x) with anticommuting coefficients. Their domain 
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on the lattice is dictated by the gauge freedom that has to be fixed. Here 
c(x),c(x) vanish at xq = L and are spatially constant at xq = 0. It is again 
convenient to extend c, 

c„(-l,p) = <v(0,p) = if (a,p)^(0,0) (4.16) 
cq(-1,0) = co(0,0), (4.17) 

and similarly for c. 

Even after gauge fixing the fluctuations, we still have covariance under 
gauge transformations of the background field. Hence the covariant differ- 
ence, like 

D^c(x) = V(x, fj,)c(x + p,)V(x, fi)' 1 - c(x) (4.18) 

in the case of c, is a natural construct to use in the following. In components 
this corresponds to 

D^x) = e^^c a (x + A) - c a (x), (4.19) 

and 

D; Ca (x) = c a (x) - e-^ x ^c a (x - p.) (4.20) 
is the corresponding covariant backward difference. 

4.3 Gluonic action 

We consider now the action (2.4) as a function of q and expand 

S(U) = S(V) + sW(q) + goS^(q) + glS^\q) + O(^ 3 ). (4.21) 

Note that the classical part <S'(V) starts at 0(gQ 2 ), but also contains higher 
order terms due to improvement, 

S(V) = [go 2 + j[c [ t 1] + g 2 ci 2) + ...]} r . (4.22) 

The simple factor 2/ L follows from the constancy of the action density of V . 
The part of the total action that is gaussian in q is given as 

S^ + S ei = (q,A iq ), (4.23) 
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where the positive linear operator Ai is defined in [3], sect. 6.4, and will 
be discussed in connection with the propagators, sect. 5 of this paper. The 
vertices of O(go,go) will now be given as a sum of terms in the form 

S w = s^ + s^ + ..., 

and analogously for other contributions. 

A particularly simple class of terms are those which only involve the 
neutral gluon components with a = 0. They are immediately generated 
from 

S(U)\ qili± =o = ~ 2 E w ^x) si* 2 [(7A- + </oW/4] (4.24) 

with the field strength /^(cc), 

f^(x) = d^q v>0 - d v q^ t o, (4.25) 

where is the forward difference operator. The relevant parts can be ar- 
ranged to give the following vertices with one to four neutral gluons: 

5 (i,a) = 2ci 1) J2M'lok/2){q k , {x)\ Xo=1 -q kfl {x)\ Xo=L _ 1 } : (4.26) 

x,fc 

S (2,a) = Ici 1 ) ^COS^/^^l^o + ^I^L-l}, (4-27) 
Z x,fc 

S {1 ' b) = ~E Ml^mfU (4-28) 

S [2 ' h) = 4E cos( 7/a ,/2)£,. (4.29) 

The contributions with q^ t ± coming from the plaquettes are numerous 
and more complicated to write down. We generated them on the computer 
(in Fortran) in a way to later allow the fast and safe contraction with the 
propagators. We here only outline the general structure and quote no lists of 
coefficients. Our example will now be the (q^,- q Vt + qx,o) vertex. For a given 
plaquette x fiu we insert into (2.4) 

tr[£/^] = tr [e 7 ^ T °e- 9o(1+D ^^ e" 90 ^ e 9oq » e 9o(1+D ^] (4.30) 
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and expand in the components q^^- The subset of terms of the desired 
structure is collected in our example in the form 

5 d,c) = £ ^e'v-.-^ + O^.+ ^ + ^l^o^ + e!). (4.31) 

The information resides in tables 4 for the directions [f r £ {/j,, u} and offsets 
£ {0,/i,z/} and in the (complex) coefflents V^ 1 ' '. Each pair of offset and 
direction is coded together into an integer which is taken as input by the 
propagator routines, when the vertex forms part of a diagram. Apart from 
the V} 1 ^ we shall also need their ^-derivatives V^'°\ which are hence stored 
in another array. The sum over r enumerates the nonzero contributions and 
has 40 terms here for plaquettes involving time and 24 terms for purely spatial 
ones without background field. Depending on the embedding of the vertex, 
the number of terms can sometimes be reduced further by the over-all reality 
of the result. 

Other vertices contributing to are constructed in the same fashion, 

5 (2,c ) = c (D E EEK (2,c V-> + Ov+> + e +)> ( 4 - 32 ) 

x =0,L-l x,fc r 

S {2 ' d) = E E Vr {2 '\,,-( X + e r)%iA X + ^ Mfli* + «£) ^,o(* + /?), 

(4.33) 

and 

(4.34) 

Each vertex has new tables although we partially use the same symbols 
here to outline the structure. It should be noted, that for terms in 
the programming efficiency is uncritical, since in 0(<7g) they will only be self- 
contracted. As we evaluate graphs in position space, the space sum is trivial 
due to translation invariance and the CPU effort for these contributions is 
negligible. 

We finally collect the term from S m here, 

^ (2 ' /) = ^E{4^,- + (^o) 2 }- (4-35) 

4 To avoid too many indices we suppress here the pLV- and Ko-dependence of these tables. 
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4.4 Ghost field action 

The Fadeev-Popov part of the action reads 

S FP = -(c,d*8 c q), (4.36) 

and 8 c q is the gauge variation of to first order in c acting as gauge parame- 
ter, d* is the covariant divergence operator, c.f. (6.13) in [3]. The expansion 
consists of the gaussian part 

S$ = (c,A c) (4.37) 
with another positive linear operator Ao and interaction vertices 

4p = - E e ' TP< W (l + \ D v) c r D^c p (4.38) 

and 

4? = i E R %,* D^D^ - q; q ; D^c a D^c T } . (4.39) 

All indices <r, r, p are summed over {0, +, — }, and e aTp is the totally antisym- 
metric symbol with e° ^ = 2i. 

5 Computation of propagators in the back- 
ground field 

In this section we discuss the definition and numerical construction of the 
propagators. This is a nontrivial problem due to the presence of a finite 
background field. 

5.1 Definition of the ghost propagator 

The ghost propagator has nonvanishing components 

G a (x,y) = {c a (x)^(y)) (no a-sum), (5.1) 
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where (.)o is the gaussian expectation value with the order zero action only. 
G a is inverse to Ao, which, with the boundary conditions quoted in the last 
section, is just the negative covariant Laplacian. This means, that 

-Y,D;D,F a (x) = f a (x) (5.2) 

holds if 

F„{x) = Y,G„{x,y)Uy) (5.3) 
y 

for any /, which fulfills the boundary conditions. G a (x,y) has to be such 
that this is also the case for F a (x). 

The symmetries of the problem suggest a spatial Fourier ansatz 

G*(*,y) = 4£e ip - (X - y) G4x°,y°;p), (5.4) 
L P 

where the p-sum runs over L 3 values suitable for the toroidal space-lattice. 
The components G a are found by solving for each p the one-dimensional 
difference equations 

|-a;a t + 6 - 2 E cos b>* + <^(*)]} ^(M'; p) = s a ,, (5.5) 

For all (<t, p) ^ (0, 0) there are equations for t,t' = 1, ... L — 1 consistent with 
the boundary conditions. This completely determines these contributions. 

The remaining component (<r, p) = (0, 0) is special, as it also contributes 
to the propagator and fulfills (5.5) at t,t' = 0, and it is fixed at this end by 
the Neumann condition (4.17). The explicit solution for this case is simply 

G o (M';0) = L -max(M')- ( 5 - 6 ) 

To compute the running coupling we need the ^-derivatives of all con- 
tributing quantities. Taking the r\ derivative of both sides in (5.5) we derive 

|-a;a t + 6 - 2 ]T cos \p k + a/?(t)]| & a (t, P ) = 
3 

-2*P'(t) sin [ Pk + *P(t)] G„(t, t'; p). (5.7) 
fc=i 

This equation defines G' a in the same way as G a as soon as the latter is 
known to furnish the right hand side instead of the ^-source in (5.5). Note 
that G' vanishes. 
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5.2 Definition of the gluon propagator 

The gluon propagator 

H^ >a (x,y) = {q fJ ,,a(x)q u a (y))o (no a-sum), (5.8) 
is inverse to Ai in (6.24), (6.25) of [3]. We directly pass to Fourier space, 

iWKy) = hA^WAy^T,^*-**'^ iW(Ay°;p), (5-9) 

L p 

with 

JO for = , u m 

^=\A W = 1,2,3 (5 - 10) 

and 

^(«) = { e ^„ for ^ =l>2>3 . I 5 ' 11 ) 

The extra phases and natural space locations at the centers of links have been 
chosen such that H^ u>(T (x ,y ;p) is real due to CP invariance (see sect. 6.3). 

The Fourier components H^^^t, t'; p) are constructed as the solutions 
Qnvit^t') (with fixed p, a) of a matrix analog of (5.5), 

{A^ K {t)Q K „(t + 1,*') + B^WQ^it, t') + A Kfl {t - l)Q w (i - 1, t')} = S^Stf. 

(5.12) 

Here and in the following summations over doubly occurring spin indices 
are implied. The coefficient matrices A,B are given in [3] by (D.4)-(D.10) 
for <T = and by (D.11)-(D.17) for a = +. They represent an invertible 
(for Ao > 0) [3] real symmetric difference operator as their components are 
real and B is a symmetric matrix. Note, that the coefficients in both cases 
depend on r\ as well as on momentum p and also on the gauge parameter 
Ao contained in S g {. The time coordinate t(t') associated with ti{y) in (5.12) 
runs over ... L — 1 for p(y)= and over 1 ... L — 1 for the other components. 
The component <r, p) = (0, 0, 0) is special in this case with the explicit 
solution 

#o,,o(M';0) = < Jo„(l + min(t,t , ))Ao 1 . (5.13) 

For the ^-derivative H 1 a one again derives an equation of the same type 
with a source term on the right hand side involving the previously constructed 
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5.3 Numerical construction of propagators 

We construct and store propagators for one pair of time coordinates (t,t') 
at a time. In a first step this is done in (spatial) momentum space. Since 
none of the three space directions is distinguished, there is symmetry under 
permuting them, which corresponds to certain discrete rotations and reflec- 
tions. Hence one gains a factor of 6 (asymptotically on large lattices) by 
solving (5.5) and (5.12) only for the reduced set of momenta pi < p2 < P3. 
Inspection of the propagator equations shows that for a = there is also 
covariance under individual reflections of any space component, which saves 
another factor of 8 in these sectors. This is not a symmetry of the full theory 
with boundary conditions, however. 

The technique of solving the difference equations for given momentum p 
will now be explained for (5.12). By two-step recursion forward and backward 
in time we construct two solutions ipj^(t) of the corresponding homogeneous 
equation 5 . They start from 

V£(-1) = V£(0) = <W (5-14) 

^L(O) = 0, ^L(l) = 4, (5.15) 

and 

tl(L) = 4,(L - 1) = <W (5.16) 

1>l(L) = 0, ^l(L-l) = 5 k „. (5.17) 

Eq.(5.12) expresses the fact that Q is the right-inverse of a certain differ- 
ence operator given by A, B. Then Q is also left-inverse, corresponding to 
the equation 

{Q^{t, t 1 + l)A VK {t') + Q^t, t')B K „{t') + Q^t, t' - l)A^{t' - 1)} = S^Stf, 

(5.18) 

whose homogeneous solutions are also given by ip^' h . If we impose now (5.12) 
and (5.18) for t ^ t' and enforce the symmetry 

Q^(t,t') = Q^(t',t), (5.19) 



5 We exclude here the case [a, p) = (0, 0) which has different boundary conditions and 
can be solved in closed form. 
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we conclude 



[tb f it)W7 l tb h At') for t < t' 

with W k a independent of t and t'. It is a consequence of our general argu- 
ments that there must exit a choice for W that makes the definition (5.20) 
consistent and solves (5.12), (5.18) for all t = t' as well. To actually determine 
W, these requirements can be shown to imply 

= 1>l(t + 1) - + 1) A XK (t) i>l(t). (5.21) 

This is a Wronskian form that is easily seen to be t-independent on the 
basis of (5.12). In the derivation one needs that ip^^t) A K \(t) ipxvft + -0 1S 
symmetric in v. This follows from its antisymmetric part being another 
Wronskian, which vanishes at the boundary due the boundary conditions. 
Also W is simplest when computed for t = which gives after some algebra 
(using (5.12) and Ajo = 0) 

Wo, = 4WL(l) + (Boo(0) + 4(-l))tt (5-22) 
Wk* = -A jk {0)^ b jv {0). (5.23) 

The computation of the ghost propagators proceeds along the same lines, 
but it is much simpler due to the absence of the spin matrix structure. 

It will be seen later, that for the most time consuming diagrams a summa- 
tion in position space will be highly advantageous (in fact mandatory to get 
to L=32). Hence, after computing and storing the momentum space prop- 
agators, we Fourier-transform them to position space after constructing all 
p-components by use of the symmetries. In position space, we then found the 
defining equations to hold essentially with machine precision. To transform 
we did not use fast Fourier transform, but the ordinary one. Therefore this 
step (for the whole set of (t, t')-values) is of computational complexity L 6 (in- 
stead of L 5 log(L) with FFT). Although the final summations are of order L 5 
at most, the propagators still consume only a small part of the CPU-time for 
L < 32. The advantage of not using FFT is the much simpler programming 
as we cannot restrict ourselves to L which are powers of two, but rather want 
all L (including the odd ones) for the final analysis of the L-dependence. In 
position space we found it more convenient to use the permutation symmetry 
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differently from before. For the gluon propagator we produce for all values of 
x — y the components fiu = 00, 01, 10, 11,12 only. The gain factor 8 for a = 
from reflections is maintained through the transformation by writing special 
Fourier routines for functions that have definite parity in all arguments. 

An initial attempt to construct propagators by (over-)relaxation directly 
in position space was abandoned. When trying to make use of the symmetries 
by identifying symmetry-related components in the iteration, we found it 
hard to still control efficiency and convergence. 

6 Summation of relevant contributions 

We now write down all contributions to V that are of 0((7q) in terms of the 
vertices to be contracted with the propagators defined in the preceding sec- 
tions. As a practical matter, we discuss issues of speed and how symmetries 
may be used to reduce the number of terms in the sum. 

6.1 Contributions to T2 

By combining now the expansion of 

5total = S + S g { + 5fP + S m (6.1) 

with (2.10) and (2.11) we have 

r J = |cP)r„ + ( 5 (L-i(s!l 1 ) 2 } o . (6.2) 

If we now exhibit the improvement structure analogously to (3.2), 

r^r^ft + Jcffr^cp^ (6.3) 

the contributions arise as follows: 

r° = (s™ + + + + 45 

-i{^) + ^) + 4 1 p ) } 2 ) , (6-4) 
c!V 2 = (S^ + S^-S^(S^ + S^ + S^)) Q , (6.5) 
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Figure 1: Expensive graphs ("Big Mac"). The dashed line represents a ghost 
propagator, the other ones gluonic fluctuations. 




Only spatial momentum zero contributes to Y c 2 , and the relevant component 
of the propagator is 

H k i,o(t, t', 0) = S kl [mm(t, t 1 ) - tt'/L]/ cos( 7ofc /2) (6.8) 

and leads to the closed expression 

T c 2 = -12sin( 7 o fc /2)tan( 7ofc /2) J L 2 ( J L - 2). (6.9) 

6.2 Efficiency and accuracy 

The computational effort to compute is by far dominated by the contri- 
butions to T% from iterated 3-point vertices that are shown in graphical form 
in fig. 1. Each of the two vertices contains a sum over the lattice. Because 
of spatial translation invariance, the total double sum has 0(L 5 ) terms when 
carried out in position space, with a prefactor of O(100 2 ) from the number 
of terms per site in each vertex. In momentum space, the complexity would 
be L 8 . A term of the structure shown in fig. 2 is 0(L 7 ) in momentum space 
and O(L) in position space which is hence clearly favoured. 

As we sum a very large number of terms we have to care about roundoff 
errors. An obvious precaution, which also leads to a more readable code, is to 



22 



Figure 2: Simple second order graphs. 



introduce a hierarchy of accumulators. Apart from this we performed the sum 
for all but the largest lattices in 32 bit accuracy in addition to the quoted 64 
bit results. The L-dependence of the discrepancy has a recognizable smooth 
trend superimposed with fluctuations, and an envelope can be estimated. We 
then expect the 64 bit results to have about 8 more significant decimal places 
than 32 bits. 

6.3 CP and T invar iance 

A possibility to still gain in efficiency is the use of further symmetries to 
reduce the number of summands. We implemented this after first verify- 
ing the symmetry relations in a full summation on smaller systems. The 
use of invariance under permutations of the three spatial directions for the 
construction of propagators was already discussed in sect. 5. 

A factor two reduction is possible by the the use of CP invariance. Action 
(2.4) and background field (2.7) are invariant under 

U{x,0) U*(x',0) (6.10) 

U(x,k) U*-\x'-k,k) (6.11) 

with x 1 = (xo, — x) for x = (so,x), which induces 6 

qoA x ) ~^ -qo-aix 1 ) (6.12) 

qkA*) -> e^^q k ,^{x'-k). (6.13) 

6 Periodicity is to be used to move x back into the range of labels chosen, xj, = . . . L—l, 
for instance. 
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and 

c a (x) -> -c_ a (x'). (6.14) 

The last relation follows from the fact that c can be viewed as a gauge 
parameter. It makes Spp CP-invariant. 

The actual reduction due to symmetries is now explained for the example 
of the left ("big mac") diagram of fig. 1. Each perturbative term in the 
position space form can be associated with its geometric origin, plaquettes 
in this case. The sum with contractions as shown can therefore be written 
schematically as 

diagram = J2 J2 /(/•*> ^ *V> u ' ^ ■> x ')- ( 6 - 15 ) 

One sees the two sets of arguments referring to the vertices. In the process 
of symmetry reduction, the various terms get different multiplicities which 
have to be incorporated in a special array and will not be mentioned in 
each case. An example is the symmetry of the diagram, which allows to 
restrict to t < t' if the non-diagonal terms are multiplied by two. Next, 
spatial translation invariance allows to pin the second plaquette to x' = 0. 
Given invariance under permutations of directions, we found it simpler not 
to reduce the x-sum, but to replace the 36 different z/, //, v' configurations 
by only 8 with appropriate multiplicities. CP corresponds to the reflection 
of the spatial coordinate of either term, which is given by the center of the 
plaquette x + (/2 + £>)/2 (c.f. (5.10)). In the sum with x' = this implies 
identical contributions from x and — x — ft — v -\- fi' -\- v' , which saves (roughly) 
a factor two in the number of different terms. 

The last symmetry to be discussed is time reversal, which is more delicate. 
We define the T-reflection as 

U(x,0) -+ t/" 1 ^' -6,0) (6.16) 
U(x,k) U(x',k) (6.17) 

where now x' = (L — xq : x). The boundaries are swapped, and the background 
field is not simply invariant. We observe however the change 

P(L-t) = -P(t)-2ir/L. (6.18) 

If we now follow the reflection by a gauge transformation 

V(x,fj,) Q(x)V(x, /j,)Q~ 1 (x + k) (6.19) 
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with 



Q, = exp (— T°(xi + X2 + xs)2ir / I?J ir\, 



(6.20) 



then V is invariant. We call this operation T. Note that Q, is antiperiodic in 
space which is admissible, as gauge fields are left periodic. 

This symmetry of the original theory is still violated by our precise way 
of gauge fixing in eq. (6.3) of [3]. As a consequence, the symmetry is not 
manifest and not immediately useful for the reduction of the sum. We hence 
replaced this condition by another possible choice, namely 

uj(x)\ X0=0 = -uj(x)\ Xo=l = k, (6.21) 

with constant k. This is immediately inherited by c, which now has to be 
constant in space with equal and opposite values at xq = 0, L. 
For the gluons we derive following [3] 

g ,o(-l, 0) = -<?o,o(£, 0) = ^(<? ,o(0, 0) - q 0fi (L -1,0)) (6.22) 

and 

#o„,o(M'; 0) = <W -\\t- Aq 1 . (6.23) 

instead of (5.13). In addition we extend c to xq = —1, L-\-l in such a way that 
D^c obeys the same boundary conditions as q^. For the ghost propagator 
(5.6) is replaced by 

G o (t,t';0) = j-^\t-t'\. (6.24) 

With these changes T invariance is manifest, and terms related by reflec- 
tions of t and t' (taking into account the natural t-location of vertices) in 
the above example become pairwise identical which (asymptotically) saves 
another factor of two. 

The stability of results under this change of boundary conditions and 
the manifest emergence of time reversal symmetry may also be viewed as an 
additional check of the computation. 
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